Search for magnetic monopoles in polar volcanic rocks 
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For a broad range of values of magnetic monopole mass and charge, the abundance of monopoles 
trapped inside the Earth would be expected to be enhanced in the mantle beneath the geomagnetic 
poles. A search for magnetic monopoles was conducted using the signature of an induced persistent 
current following the passage of igneous rock samples through a SQUID-based magnetometer. A 
total of 24.6 kg of rocks from various selected sites, among which 23.4 kg are mantle-derived rocks 
from the Arctic and Antarctic areas, was analysed. No monopoles were found and a 90% confidence 
level upper limit of 9.8 ■ 10~V 

gram is set on the monopole density in the search samples. 



The existence of magnetic monopoles was 
postulated in 1931 by Dirac as a means to ex- 
plain electric charge quantisation P, 0]. The 
Dirac quantisation argument predicts that the 
fundamental magnetic charge qm = 5ec (in this 
definition is in SI units and g is a dimension- 
less quantity) is a multiple of the Dirac charge: 
g = Ngu with g^ = 68.5 and an integer 
number. Magnetic monopoles are also funda- 
mental ingredients in grand unification theo- 
ries Although grand unification monopoles 
would typically have masses of the order of the 
unification scale (m ~ 10^^ GeV), there are 
generally no tight theoretical constraints on the 
mass of a monopole. 

Calculations within nonrelativistic quantum 
theory indicate that monopoles would bind to 
non-zero-spin nuclei through magnetic moment 
coupling, with binding energies of the order 
of several hundred kcV when assuming a hard 
core Such binding is assumed as a work- 
ing hypothesis in the present search. If iso- 
lated monopoles exist in Nature, they are sta- 
ble by virtue of magnetic charge conservation, 
and they either reside inside astronomical bod- 
ies or move freely through open space to form a 
galactic halo. Throughout this paper, "stellar" 
denotes monopoles already trapped in Stardust 
before the formation of the Solar System, and 
"cosmic" denotes free monopoles reaching the 
Solar System at a later time. 

Signatures of direct monopole pair produc- 
tion have been explored at past high-energy 
particle colliders including the LEP, HERA 
and Tevatron [sl-fioll and are being investigated 
with the Large Hadron Collider 11, l^]. How- 



ever, monopoles with masses above 7 TeV can- 
not be produced within current collider pro- 
grams. In this work, which probes monopoles 
in the mass range between the weak scale and 
the grand unification scale, it is assumed that 
monopoles may exist as relics produced out 
of thermal equilibrium in the very early Uni- 
verse. Models of cosmological inflation allow 
relic monopoles to be diluted down to noncatas- 
trophic abundances 13|. However, the vari- 



ous inflationary scenarios which have been pro- 
posed can make very different monopole abun- 
dance predictions [iJ]. Other unknowns are 
the monopole-antimonopole annihilation cross 
section and the detailed mechanisms by which 
monopoles may have bound to matter during 
primordial nucleosynthesis. Even though there 
are presently no adequate models that describe 
to which extent relic monopoles would have 
accumulated inside astronomical bodies or be 
present in cosmic rays, abundances and fluxes 
can be constrained by experiments. Monopoles 
in flight have been sought with array detec- 
tors. These set tight constraints on the flux 
of cosmic monopoles incident on Earth [H-il 
(only the most significant results are given here; 
see [29[ for a complete list). Trapped monopoles 
have previously been sought in hundreds of kilo- 
grams of samples from the Earth's crust 13 Ol- 
|35[ , in rocks from the Moon's surface 36|, |37| , 
and in meteorites 13, 2^- This work presents 
the first search for monopoles in terrestrial ig- 
neous rocks at high latitudes. 

Large planetary bodies such as the Earth 
were molten during their formation and this 
has lead to large-scale chemical differentiation. 



2 



During this early phase stehar monopoles, if 
present, will likely have sunk to the planet's 
core [38[. Stellar monopoles should therefore 
be depleted in planetary erusts, while the deep 
interiors of planets and stars, as well as the in- 
sides of some meteoroids, asteroids and comets, 
would be the only places likely to contain them 
in non-negligible amounts. 

Monopoles inside astronomical bodies of 
low viscosity possessing stable dipole magnetic 
fields would move to positions along the mag- 
netic axis where the magnetic force Fm = QmB 
{B is the vertical component of the magnetic 
field) and gravitational force Fg — ma (a is the 
gravitational acceleration) arc in equilibrium: 
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Although the early configuration of the 
Earth's internal magnetic field is poorly known, 
paleomagnetic data suggest that the Earth pos- 
sessed a dipole field since at least ~ 3.5 billion 
years The configuration of the field 



close to the Earth's core may be more com- 
plex, but the simple assumption of a dipole 
field over geologic time is reasonable. Carri- 
gan estimated that monopoles with g — go 
and m = 10^^ GeV would accumulate near the 
Earth's inner core, and developed a model of 
how monopole annihilation during geomagnetic 
reversals would contribute to the planet's in- 
ternal heat, thus limiting the grand-unification- 
mass monopole density inside the Earth to less 
than ~ 10^"'/gram [42|. On the other hand, a 
lighter mass or higher magnetic charge will raise 
the equilibrium depth. We consider monopoles 
attached to nuclei with an equilibrium position 
above the core-mantle boundary. Down to a 
depth of 2900 km, the Earth's mantle plays 
the role of an insulator between the molten 
outer core and the crust and has the proper- 
ties of a plastic solid. Although mantle dynam- 
ics are complex and various competing geody- 
namical models exist, it can generally be as- 
sumed that the mantle slowly convects as a 
whole, with a full cycle taking approximately 
400—500 million years [i^. Monopoles caught 
in the solid mantle would be unable to move 
freely. Instead, monopoles of both polarities 
would be transported up and down along with 
mantle convection regardless of the field direc- 
tion. Upon reaching the core-mantle bound- 
ary, they would sink through the liquid core 
due to the high mass, before being attracted in 
the general direction of the polar regions due 



to the magnetic charge. Over geologic time 
monopoles would migrate towards the magnetic 
axis. At the Earth's pole, a — 9.8 m-s~^ and 
B = 6.5 ■ 10~^ T, in which case Equation [1] 
yields Asurface = 1.2-10^^ GeV (presently GeV 
is a unit of mass). A monopole carrying a sin- 
gle Dirac charge [g = go) and a mass of 10^"^ 
GeV or lower would therefore be expected to 
be found beneath the Earth's polar crust and 
in melts below polar regions. A monopole car- 
rying a multiple of the Dirac charge is allowed 
to possess a proportionally higher mass. This 
mass bound is conservative because monopoles 
with equilibrium anywhere inside the mantle 
may still reach the surface through mantle con- 
vection (the core-mantle boundary corresponds 
to Aboundary = 4 • 10^^ GcV) . In a naive model, 
one may assume that monopoles would be dis- 
tributed randomly throughout the whole man- 
tle depth up to a distance from the magnetic 
axis equal to the core radius of 3400 km (this 
corresponds to latitudes > 57°), and absent ev- 
erywhere else. This results in a concentration 
of monopoles 6 times higher in polar mantle- 
derived rocks than averaged over the Earth's 
mass. 

The samples used in this search were re- 
stricted to mantle-derived igneous rocks with 
negligible levels of crustal contamination, cm- 
placed at high (> 63°) latitudes. Basaltic rocks 
from hotspots - volcanic regions under which 
the mantle is thought to be locally hotter, caus- 
ing an ascending mantle plume - are partic- 
ularly attractive as they are likely to include 
material from deep inside the mantle. Iceland 
and Hawaii are among the best known exam- 
ples of hotspots for which there is evidence 
that the erupted material comes from more 
than 600 km depth and possibly as deep as the 
core-mantle boundary [4J, |45|. Other active 
hotspot sites at high latitudes, but for which 



m- 
and 



the role of mantle plumes is debated 46 
elude Jan Mayen Island (Arctic Ocean) [47 
Ross Island (Southern Victoria Land, Antarc- 
tica) [il]. Large igneous provinces (LIPs) are 
also of interest for this work. These massive 
magmatic provinces are dominated by exten- 
sive flood basalt lavas with areal extents of 
> 100000 km^ and igneous volumes of > 100000 
km^, most of which (> 75%) was expelled dur- 
ing relatively short periods (^ 1 — 5 million 
years) (49| . Furthermore, many LIPs have been 
associated with mantle plume activity and con- 
tinental break-up [H^l. The Kap Washington 
Group volcanic sequence (North Greenland) 
and the Skaergaard intrusion (East Greenland) 
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TABLE I. Characteristics of the rock samples used in this search. If not otherwise specified, they were 
emplaced during the Cenozoic era. Control samples are indicated with (c). The latitude corresponds to 
the location at the time of emplacement. 



site 


latitude 


tectonic setting 


rock type 


samples 


mass (kg) 


Iceland [56] 


64° N 


hotspot, mid-ocean ridge 


basalt 


144 


5.916 






gabbro 


26 


1.404 


Jan Mayen Island [47] 


71° N 


hotspot 


alkali basalt 


6 


0.139 


Hawaii (c) 


21° N 


hotspot 


tholeiitic basalt 


17 


0.610 


North Greenland [5^ 


72° N 


LIP, 71-61 million 


alkali basalt, trachyte. 








years old 


trachyandesite, rhyolite 


73 


1.779 


East Greenland [58] 


68° N 


LIP, intrusion 


gabbro 


39 


1.830 


Gakljel Ridge 


84° N 


mid-ocean ridge 


tholeiitic basalt 


26 


0.707 


IVIid- Atlantic Ridge (c) 


33° S 


mid-ocean ridge 


tholeiitic basalt 


8 


0.207 


East Pacific Rise (c) 


28° S 


mid-ocean ridge 


tholeiitic basalt 


7 


0.241 


South. Victoria Land 


77° S 


hotspot 


basalt, basanite 


233 


8.163 


North. Victoria Land 


72° S 


intraplate volcanism 


basalt, trachyte 


12 


0.335 


IVIarie Byrd Land [55] 


76° S 


intraplate volcanism 


alkali basalt (HIMU) 


50 


2.184 








Iherzolite 


3 


0.148 








basalt, trachyte 


17 


0.440 


Ellsworth Land 


74° S 


intraplate volcanism 


basalt 


11 


0.300 


Horlick IVIountains 


87° S 


intraplate volcanism 


basalt 


1 


0.021 


Antarctic Peninsula (c) 


63° S 


subduction zone 


basalt 


5 


0.146 


Total search 
Total control (c) 








641 
37 


23.366 
1.204 



were considered for this search as parts of the 
High Arctic and North Atlantic LIPs, respec- 
tively 5l|, |52|. ]Vlid-occan ridges, or rift vol- 



canic zones where tectonic plates slowly move 
away from each other, are also of interest. Lava 
flows from Gakkel Ridge (Arctic Ocean) (53l.[53| 
provide attractive samples at very high lati- 
tude (84° N). Finally, some rock samples were 
selected on the basis that chemical analysis 
reveals hints of deep mantle origins. Some 
basaltic lavas from Coleman Nunatak (IVIarie 
Byrd Land, Antarctica) contain particularly 
high ^°^Pb/^"''Pb ratios (denoted as high ^, 
or HIMU), which indicates low extent of melt- 
ing and relatively deep origin [H^ . In addition, 
some of the lavas carry nodules of Iherzolite, 



which have been carried up from the mantle 
source rocks without melting. Control samples, 
which should not contain stellar monopoles be- 
cause they fail one of the search criteria, were 
also included: crust-derived lavas from a sub- 
duction zone (Antarctic Peninsula), and sam- 
ples from a hotspot or mid-ocean ridge at low 
latitude (Hawaii, IVIid- Atlantic Ridge and East 
Pacific Rise). The samples were shaped ei- 
ther as cylinders of 2.5 cm diameter and about 
2.5 cm length, or crushed into fragments, which 
were placed into plastic cuboid boxes 2.3 cm on 
one side. The analysed samples are listed in Ta- 
ble U and amount to a total of 23.4 kg of search 
samples and 1.2 kg of control samples. 



Samples were measured with a 2G Enter- 
prises, model 755R, 3-axis DC-SQUID rock 
magnetometer housed in a shielded room at 
the Laboratory of Natural Magnetism, ETH 
Zurich. For magnetic dipoles the current re- 
verts to zero on complete passage through the 
magnetometer superconducting coils. However, 
a monopole would leave the signature of a per- 
sistent current. This technique allows us to di- 
rectly measure the magnetic charge contained 
inside a sample without the need to extract 
monopoles and with no mass dependence. Cur- 
rent measurements were performed in steps, in- 



cluding measurements where the sample is in- 
side the sensing coils as well as 50 cm away from 
the sensing coils before and after the pass. Oc- 
casional passes with an empty sample holder 
were made for background subtraction. The 
persistent current is defined as the measured 
value after pass minus the value before pass 
(subtracting the same quantity for the empty 
holder) , normalised such as to give the strength 
of magnetic pole contained in the sample in 
units of gn- As described in detail in [59|, cal- 
ibration was performed using the convolution 
method, which consists of profiling the magne- 
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tometer response as a function of distance for a 
sample with well-known magnetisation and in- 
ferring the response for a monopole. As a cali- 
bration cross-check, the response to a magnetic 
pole was tested by introducing one extremity 
of a thin solenoid of 25 cm length with applied 
currents corresponding to values of magnetic 
charge of 0.124 go, 1.24 go, 12.4 go and 124 go- 
The two methods yield consistent results within 
a normalisation uncertainty of 10%. 

Samples with a total magnetisation > 1.5 • 
10^ go (or magnetic dipole moment > 4.4 ■ 
10~^ Am^) were found to sometimes cause the 
flux-lockcd loop of the SQUID to be lost and 
recovered at a different quantum level. This 
leaves a signal similar to what is expected from 
a monopole. Weaker moments generally did 
not show this effect. Precautions were there- 
fore taken so that all samples would have mag- 
netisation levels below 1.5 ■ 10^ go- Crushing 
the sample material into a gravel- or sand-sized 
powder randomises the magnetic moments from 
the constituent ferromagnetic minerals, which 
reduces the dipole signal. This method was fre- 
quently used in this study. Alternatively, the 
magnetisation can be reduced by more than an 
order of magnitude by exposing the sample to 
an alternating field. There is no risk of dislodg- 
ing a trapped monopole if a binding energy of 
100 kcV or more is assumed. Demagnetisation 
was carried out only on 10% of the Antarctic 
samples probed in this study. 

Measurements of persistent currents after 
first passage through the magnetometer are 
shown for all samples in Fig. [T] (top). In the 
range from —0.1 to 0.1 go, the distribution is 
Gaussian with mean value —0.002 ± 0.001 go 
and standard deviation 0.026 ±0.001 go- Non- 
Gaussian tails slightly extend the distribution 
beyond this range. Five candidates out of 
678 samples yield absolute values which devi- 
ate from zero by more than 0.25 go- The two 
first of these candidates yield the largest val- 
ues (0.8 go and 1.6 go) and also have total 
magnetisations in excess of 10^ go, close to the 
1.5 • 10^ go limit beyond which measurements 
are known to be unreliable. Additional mea- 
surements of the five candidates using various 
orientations of the samples arc shown in Fig. [1] 
(bottom). These multiple measurements con- 
firm the zero magnetic charge hypothesis. It is 
possible to get a rough estimate of the proba- 
bility that a random sample containing a gen- 
uine monopole with \g\ = gu would yield a per- 
sistent current close enough to zero to remain 
unnoticed. The probability to mismeasure the 
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FIG. 1. Top: persistent current after first passage 
through the magnetometer for all samples. Bot- 
tom: results of repeated measurements of candidate 
samples with absolute measured values in excess of 
0.25 go- 



current by an absolute value which deviates 
from go by less than 0.25 go is about 0.3% 
(out of 678 samples, only the first candidate dis- 
cussed above satisfies this condition, but some 
of the other candidates are close enough that we 
conservatively assume two). The probability to 
mismeasure the current in the direction where 
it would cancel out the current induced by a hy- 
pothetical monopole (whose charge can be pos- 
itive or negative) is 1/2. Thus we obtain that 
0.3%/2 = 0.15% of the signals with = go 
would escape detection; less ii \g\ > go- It 
is concluded that no monopoles with magnetic 
charge \g\ > go were present in the samples. 

The most extensive meteorite search to date 
- the only other direct search with a non- 
negligible sensitivity to stellar monopoles - sets 
a limit on the monopole density in meteoritic 
material of less than 2.1 • 10^^/gram at 90% 
confidence level. The study analysed 112 kg 
of meteorites [sHj, among which 100 kg are 
chondrites and can thus be assumed to consist 
of undifferentiated material from the primary 
solar nebula. This represents a little more than 
4 times more material than used in the present 
search. As discussed above, for monopole mass 
and charge satisfying Equation [1] for a posi- 
tion above the core-mantle boundary, this dif- 
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ference can be compensated for by an increase 
in monopolc concentration of roughly a factor 6 
in polar mantle-derived rocks due to monopole 
accumulation along the Earth's magnetic axis. 
One can think of two ways in which these re- 
sults on stellar monopoles could be further im- 
proved in the future: by probing large (> 100 
kg) amounts of meteorites and polar rocks with 
a high-efficiency magnetometer, or by gaining 
access to new types of samples such as asteroid 
and comet fragments. 

In summary, massive monopoles of stellar ori- 
gins would be absent from planetary surfaces 
and would tend to accumulate along the mag- 
netic axis in planets with internal magnetic 
fields. If monopoles in the mass range 10"^ < 
m < 10^"^ GeV are present within the Earth, 
they would be expected to be found inside the 
Earth's mantle below the geomagnetic poles. 
Assuming that monopoles bind strongly to nu- 
clei, they would be trapped in mantle-derived 
rocks. This paper presents the first search for 
monopoles in polar igneous rocks. The search 
probed 23.4 kg of samples, for which a limit 
on the monopole density of 9.8 ■ 10~^/gram at 
90% confidence level is set, which in a simple 
model translates into a limit of 1.6 • 10~^/gram 
in the matter averaged over the whole Earth. 
This search has a comparable or better sensi- 
tivity than the most extensive meteorite search 
and provides a novel probe of stellar monopoles 
in the Solar System. 
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the United States Polar Rock Repository, which 
is sponsored by the United States National 
Science Foundation, Office of Polar Programs. 
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the Swiss National Science Foundation and a 
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